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Abstract 

We consider the homogenization of the Navier-Stokes equation, set in a channel with 
a rough boundary, of small amplitude and wavelength e. It was shown recently that, 
for any non-degenerate roughness pattern, and for any reasonable condition imposed at 
the rough boundary, the homogenized boundary condition in the limit e = is always 
no-slip. We give in this paper error estimates for this homogenized no-slip condition, and 
provide a more accurate effective boundary condition, of Navier type. Our result extends 
those obtained in [5J [T3] , in which the special case of a Dirichlet condition at the rough 
boundary was examined. 
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1 Introduction 

Most works on Newtonian liquids assume the validity of the no-slip boundary condition: the 
velocity field of the liquid at a solid surface equals the velocity field of the surface itself. 
This assumption relies on both theoretical and experimental studies, carried over more than 
a century. 

Still, with the recent surge of activity around microfluidics, the question of fluid-solid 
interaction has been reconsidered, and the consensus around the no-slip condition has been 
questioned. Several experimentalists, observing for instance water over mica, have reported 
significant slip. More generally, it has been claimed that, in many cases, the liquid velocity 
field u obeys a Navier condition at the solid boundary S: 

{I d - v® v)u\x = X(I d - v® v)D{u)v\^ y u-i/|e = 0, A>0 (Na) 

where v is an inward normal vector to S, and D{u) is the symmetric part of the gradient. Slip 
lengths A up to a few micrometers have been measured. This is far more than the molecular 
scale, and would therefore invalidate the (macroscopic) no-slip condition 

«|e = 0. (Di) 

Nevertheless, such experimental results are widely debated. For similar experimental settings, 
there are huge discrepancies between the measured values of A. We refer to the article [18J 
for an overview. 
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In this debate around boundary conditions, the irregularity of the solid surface is a major 
issue. Again, its effect is a topic of intense discussion. On one hand, some people argue 
that it increases the surface of friction, and may cause a decrease of the slip. On the other 
hand, it may generate small scale phenomena favourable to slip. For instance, some rough 
hydrophobic surfaces seem more slippery due to the trapping of air bubbles in the humps of 
the roughness. Moreover, irregularity creates a boundary layer in its vicinity, meaning high 



velocity gradients. Thus, even though (Di) is satisfied at the rough boundary, there may be 
significant velocities right above. In other words, the no-slip condition may hold at the small 
scale of the boundary layer but not at the large scale of the mean flow. This phenomenon, 
due to scale separation, is called apparent slip in the physics litterature. 

In parallel to experimental works, several theoretical studies have been carried, so as 
to clarify the role of roughness. Many of them relate to homogenization theory. First, 
the irregularity is modeled by small-scale variations of the boundary. Then, an asymptotic 
analysis is performed, as the small scales go to zero. The idea is to replace the constitutive 
boundary condition at the rough surface by a homogenized or effective boundary condition at 
the smoothened surface. In this way, one can describe the averaged effect of the roughness. 
We stress that such homogenized conditions (often called wall laws) are also of practical 
interest in numerical codes. They allow to filter out the small scales of the boundary, which 
have a high computational cost. 

Let us recall briefly the main mathematical results on wall laws. To give a unified descrip- 
tion, we take a single model. Namely, we consider a two-dimensional rough channel 

n £ := QUY,UR £ 



where O = Mx(0,l) is the smooth part, R £ is the rough part, and S = Kx {0} their interface. 
We assume that the rough part has typical size e, that is 

R* : = £ R, R := {y, > y 2 > oj { Vl )} 

for a Lipschitz function u : R i— > (— 1, 0). We also introduce 

r £ := e r, r := {y, y 2 = u fa)} 

See Figure [T] for notations. We consider in this channel a steady flow u £ . It is modeled by the 
stationary Navier-Stokes system, with a prescribed flux <p across a vertical cross-section a 6 of 
Q £ . Moreover, to cover all interesting cases, we shall consider either pure slip, partial slip or 
no-slip at the rough boundary T e . This means that the constant A below shall be either +oo, 
positive or zero. For simplicity, we assume no-slip at the upper boundary. We get eventually 

-Au £ + u £ ■ Vu £ + Vp £ = 0, x G n £ , 
div u £ = 0, x £ Q e , 



< 




(NS E ) 



„ (Id — v u)u £ \r^ = \ £ (Id— v ® v)D(u £ )v\ri:, u £ • = 0. 

Notice that the flux integral in the third equation does not depend on the location of the cross- 
section o~ £ , thanks to the divergence- free and impermeability conditions. We also emphasize 
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Figure 1: The rough domain Q e . 



that this problem has a singularity in e, due to the high frequency oscillation of the boundary. 
Thus, the problem is to replace the singular problem in fl £ by a regular problem in Q. The 
idea is to keep the same Navier-Stokes equations 



-Aw + u ■ Vu + Vp = 0, x G ft, 
div u = 0, x G fl, 



u 



X2 = l 



0. 



(NS) 



but with a boundary condition at the artificial boundary S which is regular in e. The problem 
is to find the most accurate condition. 

A series of papers has addressed this problem, starting from the standard Dirichlet con- 
dition at T £ (A e = in (NS £ )). Losely, two main facts have been established: 



1. For any roughness profile oo, the Dirichlet condition (Di) provides a 0(e) approximation 

°f u£ in L lioM- 

2. For generic roughness profile u, the Navier condition does better, choosing A = ae for 



some good constant a in (Na). 



Of course, such statements are only the crude translations of cumulative rigorous results. 
Up to our knowledge, the pioneering results on wall laws are due to Achdou, Pironneau and 
Valentin [U [3], and Jager and Mikelic [141 I15j . who considered periodic roughness profiles 
oo. See also p£] on this periodic case. The extension to arbitrary roughness profiles has been 
studied by the second author (and coauthors) in articles [6j [T2], [13] . The expression generic 
roughness profile means functions oo with ergodicity properties (for instance, oo is random 
stationary, or almost periodic). We refer to the forementioned works for all details and 
rigorous statements. Let us just mention that the slip length ae is related to a boundary 
layer of amplitude e near the rough boundary. It is the mathematical expression of the 
apparent slip discussed earlier. 

Beyond the special case A e = 0, some studies have dealt with the general case A e G [0, +oo]. 
The limit u° of u 6 , and the condition that it satisfies at £ have been investigated. In brief, 
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the striking conclusion of these studies is that, as soon as the boundary is genuinely rough, 
u° satisfies a no-slip condition at E. This idea has been developped in [9] for a periodic 
roughness pattern uj. It has been generalized to arbitrary roughness pattern in [8]. In this 
last article, the assumption of genuine roughness is expressed in terms of Young measure. 
When recast in our 2D setting, it reads: 

(H) The family of Young measures (dn yi ) associated with the sequence (u/(-/e)) g is s.t. 
dfi yi ^ 5o (the Dirac mass at zero), for almost every y\ E R. 
Under (H), one can show that u £ locally converges in ff 1 -weak to the famous Poiseuille flow: 
u°(x) = (U°(x 2 ),Q), U°(x 2 ) = 6<j)x 2 (l-x 2 ) 



which is solution of (NS )-(|Di|). We refer to [8j for all details. 



This result can be seen as a mathematical justification of the no-slip condition. Indeed, 
any realistic boundary is rough. If one is only interested in scales greater than the scale e 



of the roughness, then (Di) is an appropriate boundary condition, whatever the microscopic 
phenomena behind. Still, as in the case A e = 0, one may be interested in more quantitative 
estimates. How good is the boundary condition? Can it be improved? Is there possibility 
of a 0(e) slip? Such questions are especially important in microfluidics, a domain in which 
minimizing wall friction is crucial (see |23j). 

The aim of the present article is to address these questions. We shall extend to an arbitrary 
slip length A £ the kind of results obtained for X £ = 0. Of course, as in the works mentioned 
above, we must assume some non-degeneracy of the roughness pattern. We make the following 
assumption: 

(H') There exists C > 0, such that for all 2-D fields u € lR) satisfying u • u\r = 0, 

IMIl 2 (R) < C \\Vu\\ L 2( R y 

Assumption (H'), and its relation to the assumption (H) will be discussed thoroughly in the 
next section. Broadly, we obtain two main results. The first one is 



Theorem 1. There exists 4>$ > 0, such that for all \(j>\ < <fio, for all e < 1, system (NS £ ) has 
a unique solution u e in Hl loc (tt £ ). Moreover, if \ £ = or if (W) holds, one has 



u 



"'Iff 1 , (n«) ^ C4>y/e, ||-u £ -u°|| L 2 (n) < C<f)e, 



where u is the Poiseuille flow, satisfying (NS)-(Di). 



In short, the Dirichlet wall law provides a 0((pe) approximation of the exact solution u £ in 
^mZoc(^)' f° r an y ^ 6 ^ [0, +oo]. This gives a quantitative estimate of the convergence results 
obtained in the former papers. Note that the dependence of the error estimates on both <fi 
and e is specified. In the case X £ = 0, this improves slightly the result of [6], where the <f> 
dependence was neglected. 

Our second result is the existence of a better homogenized condition. Here, as outlined 
in article [13], some ergodicity property of the rugosity is needed. We shall assume that uj is 
a random stationary process. Moreover, we shall need a slight reinforcement of (H'), namely: 
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(H") There exists C > 0, such that for all 2-D fields u £ (i?) satisfying u ■ u\y = 0, 

\\u\\ L 2 {R) < C\\D(u)\\ L 2 {R) , D(u) = -(Vu + (Vu) t ). 
We shall discuss this assumption in section [2] We state 

Theorem 2. Let uj be an ergodic stationary random process, with values in (—1,0) and K- 
Lipschitz almost surely, for some K > 0. Assume either that X s = 0, or that X s = A > for 
all e, and the non- degeneracy condition (H") holds almost surely, with a uniform C. Then 



there exists a > and <j>Q > such that, for all < 4>q, e < 1, the solution u N of (NS)-(Na 
with X = a e satisfies 



\ 1/2 

sup — / \u e — u N \ 2 dx\ = o(e), almost surely. 



1 

k r>1 R 7nn{|x-i|<i?} 

We quote that the norm above is common in the framework of stochastic pde's: see for 
instance [5]. We also quote that, even in the case X s = 0, this almost sure estimate is new: the 
estimates of [6] involved expectations. This result can also be extended to other slip lengths 



A e in (NS £ ); more precisely, up to a few minor modifications, our techniques also allow us to 
treat slip lengths A £ such that A £ 3> 1, or A e < e 2 , or A e = A°e. 

Briefly, the outline of the paper is as follows. In section [2j we will discuss in details the 
hypotheses (H') and (H"). Section [3] will be devoted to the proof of theorem [TJ In section |4j 
we will analyze the boundary layer near the rough boundary. This will allow for the proof of 
Theorem [2j to be achieved in section [5} 



2 The non-degeneracy assumption 

The goal of this section is to discuss hypotheses (H') and (H"), and, in particular, to give 
sufficient conditions on the function u for (H') and (H") to hold. We will also discuss the 
optimality of these conditions in the periodic, quasi-periodic and stationary ergodic settings, 
and compare them to assumption (H). 

2.1 Poincare inequalities for rough domains: assumption (H') 

First, let us recall that if the non-penetration condition u • v\r = is replaced by a no-slip 
condition u\r = 0, then the Poincare inequality holds: indeed, for all u £ H l (R) such that 
u\y = 0, we have 



\u{y\,m)\ 2 dy\ dy 2 



R 



R 



Hi 



d 2 u(y 1 ,t)dt 



dyi dy 2 



< C I \d 2 u{y 1 ,t)\'dy 1 dt, 

IR 



where the constant C depends only on ||a;||x,oo . 

Assumption (FT) requires that the same inequality holds under the mere non-penetration 
condition; of course, such an inequality is false in general (we give a counter-example below in 
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the case of a fiat bottom). In fact, (H') is strongly related to the roughness of the boundary: if 
the function uj is not constant, then the inward normal vector v = (1 + w /2 ) _1 / 2 (— u/, 1) takes 
different values. Since u ■ v\r = 0, we have a control of u in several directions at the boundary 
(at different points of T). In fine, this allows us to prove that the Poincare inequality holds, 
and the arguments are in fact close to the calculations of the Dirichlet case recalled above. 
• We now derive a sufficient condition for (H'): 

Lemma 3. Let co 6 W 1,co ()R.) with values in (—1,0) and such that supw < 0. Assume that 



(■A 

3,4 >0, inf / b'fyi + t)\ 2 dt > 0. 
2/ieR./o 



(2.1) 



Then assumption (H') is satisfied. 
Proof. The idea is to prove that for some well-chosen number B > 0, there holds 

[ \u(y)\ 2 dy < C B [ [ \u(y 1 ,y 2 )-u(y 1 +t)\ 2 dtd Vl dy 2 (2.2) 
JR JrJo 

< C B [ \Vu(y)\ 2 dy. (2.3) 
JR 

The first inequality is a direct consequence of assumption (2.1). The proof of the second one 
follows arguments from [9] , and is in fact close to the proof of the Poincare inequality in the 
Dirichlet case. 

First, for all B > 0, we have 

\u(yi,y2) • v{yi + t)\ 2 dt dy 1 dy 2 



R Jo 




B 



1 



> 



1 + W 



(-ui(y)u' (yi + t) + u 2 {y)) dt dy\ dy 2 

B 



R J l + u' 2 ( yi +t) 

[ dyu\{y) [ J{ yi + t) 2 dt + B [ u 2 

JR JO JR 



> 



ui{y)u 2 (y){uj{yi + B) - co(yi)) dy 

1 inf ( B - IHloo, inf f (J(y\ + t) 2 dt - \\u 
1 + IMI« V vi^Jo 



\u(y)\ 2 dy. 



R 



Assume that B > A, and set 



a := inf 

Jo 



\u'( yi + t)\ 2 dt. 



Notice that a > thanks to (2.1). Then 

"B 



inf 

2/1 6M Jo 



J( yi + tydt > 



B 
A 



a, 



and thus there exists a positive constant c such that for all B > A, 

/ / Kyi, 2/2) • v(y x +t)\ 2 dtdyi dy 2 > c(B - 1) / \u(y)\ 2 dy. 
JRJO JR 
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Thus for B large enough, inequality (2.2) is satisfied. 

As for (2.3), let us now prove that for all B > 0, there exists a constant Cb such that 

n\u{yi,y2)-v{yi + t)\ 2 dtdyidy 2 <C B / \Vu(y)\ 2 dy. 
JR 

We use the same kind of calculations as in [9]. The idea is the following: for all y £ R, 
t G [0,B], let 

z = (y 1 +t,uj{y 1 + t))eT. 

Let £y )t be a path in W l >°°{[0, 1],K 2 ) such that l yjt (0) = y, £ y>t (l) = z and ty )t {r) G for all 
r € (0,1). Then 

u(y)-u(z) = [ « )t (r).V)«MT))dT, 

JO 

and thus, since u(z) • z/(yi + i) = 0, 

|n(y)-^( yi + t)| < / |(C(r)- V)n(^(r))| drdtdy. 
Jo 

There remains to choose a particular path lyj. 

Notice that in general, we cannot choose for the straight line joining y and z, since 
the latter may cross the boundary T. We thus make the following choice: for A G (supw,0), 
we set 

A ■= (yi.A), 

^ := (yi+i,A). 

We define the path £ V)t by 

^,*(0) = y, £ y ,t = 4. V ( f ) = Z A, V (!) = 2, 

and £ yi t is a straight line on each segment [0, 1/3], [1/3, 2/3], [2/3, 1] (see Figure [2]). 

Notice that £ Vj t depends in fact on A, although the dependance is omitted in order not to 
burden the notation. With this choice, we have 



\u(y) ■ v{yi + t)\ < / \d 2 u\{yi,y' 2 ) dy' 2 

■%2,A] 



< 



+ I \d 1 u\(y 1 +y' 1 ,\)dy' 1 
Jo 

+ / \9%u\{yt + t,y' 2 ) dy 2 

r-B 

I3H (2/1,2/2) d y'2 + I \d\u\{yi + y[A) dy[ 

0,(2/1) Jo 





+ / \<hv\(yi + t,y' 2 ) dy 2 . 

J ui(vi +t) 



lui(y 1 +t) 

Integrating with respect to y and t, we obtain, for all A G (sup u, 0) 

|2 



\u(y)-v( yi +t)\ A dydt<C B [ / \d 2 u\\y)dy+ / \d\u\ {y\, A) dyi 

y&RJO \JR JR 

Integrating once again with respect to A yields the desired inequality. □ 
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j-0 



t.'=(Y\ .lambda) 



z " = (yi +t ' lambda) 




z-(yi+t. omega(7!+t)) 



Figure 2: The path £ y j. 



Let us now examine in which case assumption (2.1) is satisfied in the periodic, quasi 



periodic and stationary ergodic settings: first, if uj is T-periodic, where T := M/Z, then (2.1) 
merely amounts to 



J 2 > 0. 



Hence (H') holds as soon as the lower boundary is not flat. In this case assumption (2.1) is 
necessary, as shows the following example: assume that u = — 1, and consider the sequence 
(uk)k>i m ^(R) defined by u^^ = and 



Uk,i(y) 



l if|yi|<fc, 
if|yi|>fc + l, 



and Uk,i(y) G [0, 1] for all y, \\Vuk,i\\L°° < 2. 

Then it is easily checked that Uf. ■ v\t = 0, and that 



\ u k\\L 2 (R) 



> 2k. 



On the other hand, 



IIV^H 2 



L 2 (R) 



\Vu k \ 2 < 8 V/s > 1. 

'k<yx<k+l 

Hence assumption (H') cannot hold in R. 

In the quasi-periodic case, the situation is similar to the one of the periodic case, i.e. 



(2.1) ^=> J ^ 0. 



Indeed, assume that 

uj{y x ) = F(Xyi) 
for some A G M. d , F £ C 2 (T d ), with d > 2 arbitrary. Then 

[ A oj' 2 ( yi +t)dt= [ A (X ■ VF) 2 (X( yi + t)) dt. 
Jo Jo 
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Write F Fourier series: 



F(Y) = ^ a k e 2ink - Y VY G T d . 



Then 



rA 



(A • VF) 2 (X(yi + t)) dt = -2vr a k ai (A • fc) (A • l)e 2i 



in(k+l)-Xy 



k,l& 
A-(fc+Z)/0 



+4vr 2 A a k ai(X-k) 2 . 



k,l£Z d , 
A-(fc+/)=0 

The first term is bounded uniformly in y\ and A provided the sequence a k is sufficiently 
convergent and A satisfies a diophantine condition. Consequently setting 



C = 4vr 2 a k ai(X-k) 2 , 



k,iez d , 
\-(k+l)=0 

we deduce that there exists a constant C such that 

-A 



r 

MA > 0, Vyi G M, C A - C < oo' 2 (yi + t) dt < C A + C. 

Jo 

The above inequality entails that Co > 0. If Co > 0, inequality (2.1) is proved. If Co = 0, we 
infer that 



' 2 ^ 
u < oo. 



As a consequence, since u/ is uniformly continuous on M, lim^i^oo w'(t) = 0. On the other 
hand, it can be proved thanks to classical arguments that for all e > 0, N > 0, there exists 
n G N such that n > N and 

d(\n,Z d ) < e. 

For e small and N large, and t in a fixed and arbitrary bounded set, we obtain 

u)'(t + n) = o(l) 

= A • \7F(Xt + An) 

= A • VF(Xt) + o(l) 

= c/(t) + o(l). 

Thus w'(t) = 0, and w' = 0. 



Hence we deduce that (2.1) is satisfied as soon as u' is not identically zero, at least for 
"generic" quasi-periodic functions (i.e. such that the Fourier coefficients of the underlying 
periodic function are sufficiently convergent and such that A satisfies a diophantine condition) . 
In fact, slightly more refined arguments (which we leave to the reader) show that the result 
remains true as long as 

X! \ k \ \ a k\ < °°> 

fcez d 
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without any assumption on A. 

Let us now give give another formulation of (2.1) in the stationary ergodic case. We 
denote by (M, //) the underlying probability space, and by (T yi ) yi £R the measure-preserving 
transformation group acting on M. We recall that there exists a function F G L°°(M) such 
that 

w(yi,m) = F(r yi m), y\ G R, me M. 
As in [6] , we define the stochastic derivative of F by 

d m F(m) := u/(0, to) Vto G M, 

so that u)'(yi,m) = d m F{T yi m) for (yi, m) elx M. We claim that almost surely in to G M, 

/■A /-A 

inf / |u/(yi + t,m)| 2 = essinf m / e a/ / |<9 m F(TiTO')| 2 dt. 



.yi 



Indeed, notice that the left-hand side is invariant under the transformation group (t z1 ) Z1 £r 
as a function of to G M. As a consequence, it is constant (almost surely) over M; we denote 
by (ft the value of the constant. Since uj' G L°°, we also have 



yi 

Now, for all y\ G Q, 



f A 

inf / |w'(yi + t, m)| 2 (it a.s. in M. 

'ieQ Jo 



w'(yi + t,TO,)| 2 dt = / \d m F(T t (T yi m)) dt 



nA 

> essinf m / eA f / \d m F(T t m')\ 2 dt =: eft' 
Jo 



almost surely in M. Taking the infimum over y\ G Q, we infer that (ft > 0'. 

On the other hand, by definition of (ft' , for all e > 0, there exists M £ C M such that 
P{M £ ) > and 



^ < / \d m F(T t m)\ 2 dt<(ft' + e VmeM f 
Jo 



Consequently, for all m G .M e , we have 

».4 



that is, 



/■A /-A 

inf / \u'(yi +t,m)\ 2 dt < / \<J(t, m)\ 2 dt < (ft' + s, 
nmJo Jo 



<<ft' + e. 



Hence (ft = (ft' . Eventually, we deduce that in the stationary ergodic case, assumption (2.1) is 
equivalent to 

f A 

3A > 0, essinf m£ A/ / \d m F(T t m)\ 2 dt > 0. (2.4) 
Jo 

A straightforward application of the stationary ergodic theorem shows that (2.4) implies that 

E\\d m F\ 2 } > 0. 
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However, assumption (2.4) appears to be much more stringent than the latter condition: 
indeed, (2.4) is a uniform condition over the probability space M, whereas the convergence 



1 f R 

- / \d m F(r t m)\ 2 dt — > E[\d m F\ 2 } 



only holds pointwise. 

• Let us now compare condition (2.1) with the assumption (H) of [8]. In order to have 
a common ground for the comparison, we assume that the setting is stationary ergodic. In 
this case, the family of Young measures associated with the sequence u/(-)/e can be easily 
identified: indeed, according to the results of Bourgeat, Mikelic and Wright (see [7]), for all 
G £ C X (K) and for all test function <p E L X (R x M), there holds 



/ 



G (u)' ( — ,mj ) <p(yi,m) dyi d/i(m) -> / E[G(d m F)](f(y 1 ,m) dyi dfi(m). 

xM v v e 77 JlxM 



By definition of the Young measure, the left-hand side also converges (up to a subsequence) 
towards 

/ (G,dfi yi )(p(yi,Tn))dyidfj,(m). 

JRxM 

As a consequence, we obtain 

(G, dfi yi ) = E[G(d m F)\ for a.e. y x £ M. 
Hence condition (H) is equivalent (in the stationary ergodic setting) to 

-E[|<9 m -F| 2 ] > 0, i.e. F non constant a.s. 

We deduce that assumptions (H) and (H') are equivalent in the periodic and quasi-periodic 
settings. In the general stationary ergodic setting, however, condition (2.1) is stronger than 
(H). But since we do not know whether (2.4) is a necessary condition for (H') in the stationary 
setting, we cannot really assert that (H') is stronger than (H). 



2.2 Korn-type inequalities: assumption (H") 

We now give a sufficient condition for (H"). Notice that our work in this regard is related to 
the paper by Desvillettes and Villani [10], in which the authors prove that for all bounded 
domains $7 C M. N which lack an axis of symmetry, there exists a constant K{Q) > such that 

\\D{u)\\ L 2 m > K (fi)||Vu||^ ( n) 
Vu G H l {Q) N s.t. u-v\ dn = 0. 

The differences with our work are two-fold: first, in our case, the domain R is an unbounded 
strip, which prevents us from using Rellich compactness results in order to prove (H"). More- 
over, the tangency condition only holds on the lower boundary of R. However, as in |10| . 
we show that condition (H") is in fact related to the absence of rotational invariance of 
the boundary T. Let us stress that this notion is related, although not equivalent, to the 
non-degeneracy assumption of the previous paragraph (see ( |2.1[ )). 
We first define the set of rotational invariant curves: 
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Definition 4. For (yo,R) E M 2 x [0, oo), denote by C(yo,R) the circle with center yo an d 
radius R. 

For all A > 0, we set 

K A := { 7 E ^^([0,i]) 2 , 3(1/0, i?) G M 2 x (0,oo), 7 ([0,A]) cC(y ,#)} 
U{ 7 E W rl ' DO ([0,A]) 2 ,i/ 7 = est}, 

where v~ is a normal vector to the curve 7, namely 

1 /W 2 



7 (7i 2 + 7 ^)V7i 

Notice that 7£a is a closed set with respect to the weak - * topology in W 
We then have the following result: 

Lemma 5. Let co E W 1,oc (R). For A> 0, k E Z, let 

7^ : yi E [0,A] ^ (y u u)(y x + kA)). 
Assume that there exists A > such that 



{7jf,fc€Z}n^A = ) (2-5) 
where the closure is taken with respect to the weak - * topology in W ,ao . Then (H") holds. 



Assumption (2.5) means that each slice of length A of the boundary remains bounded 
away from the set of curves which are invariant by rotation. In particular, in the periodic 
case, a simple convexity argument shows that all non-flat boundaries satisfy (2.5) (it suffices 
to choose A equal to the period of the function oS). 

The proof of Lemma [5] uses the following technical result: 

Lemma 6. For all Y > sup a;, let 

R Y := {y E M 2 , u(y x ) < y 2 < Y}. 

Consider the assertion 

(K Y ) 3Cy > 0, Vu E H 1 (Ry) s.t. u ■ v\y = 0, j \u{y)\ 2 dy < C Y ! \D{u)\ 2 . 

J Ry ^ Ry 

If there exists Yq such that (Ky ) is true, then (Ky) is true for all Y > sup a;. 
We postpone the proof of Lemma [6] until the end of the section. 

Let us now prove Lemma [5] the idea is to reduce the problem to the study of a Korn-like 
inequality in a fixed compact set, and then to use standard techniques similar to the proof of 
the Poincare inequality in a bounded domain. 

First step: reduction to a flat strip. 

According to Lemma [HJ it is sufficient to prove the result in a domain Ry for some 
Y > supcj sufficiently large (notice that the boundary T is common to all domains Ry). 

We use the extension operator for Lipschitz domains defined by Nitsche in |20| . Since the 
result of Nitsche is set in a half-space over a Lipschitz curve, we recall the main ideas of the 
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construction, and show that all arguments remain valid in the case of a strip, provided the 
width of the strip is large enough. 

We denote by £1- the lower half-plane below T, namely 

:= {y G R 2 , y 2 < wfcft)}. 

According to the results of Stein (see [22J ) , there exists a "generalized distance" 6 G 
C°°(0_) such that 

< 2(w(yi) - y 2 ) < S(y) < C (oj( yi ) - y 2 ) Vy G fi_, 
\d%8{y)\ < C^iy) 1 -^ Va G N 2 Vy G 

(In general, since cj is merely a Lipschitz function, the function d(-,T) has very little regularity, 
whence the need for a generalized distance.) 
Let ip G C([l, 2]) such that 

^ ^(A) dA = 1, ^ AV»(A) dA = 0. 

For u G H 1 (Ry), define an extension u of u in a strip (info; — r/, y) for some r\ > by 
w(y) = w(y) ifyGiZy, 

: = J V'(A) [^(z/a) + Ad i( 5(y) « 2 (yA)] dA, if y G 0_ 
where y A := (yi,y 2 + A5(y)), 



Choose Y such that 

F > 2C supw - (2C - l)infw. 

Then if rj > is sufficiently small, y\ G Ry for all y G (infw — r/, Y). The function u thus 
defined does not have any jump across T. Moreover, it can be checked that 

l|£(«)IU2(Rx(iiifu;-J7,y)) < C \\ D (u)\\ L 2( R y 
Indeed, if y G 0_ and y 2 > inf u — r], 

[diUj + dj-Ui] (y) = J dX^(X) [(cVitj + (y A ) + 2X 2 d i 5(y)d j 5(y)d 2 u 2 (y x ) 

+Xdi5(y) {d 2 Uj + dju 2 ) (y A ) 
+A<9,£(y) (9 2 ^i + diu 2 ) (y A ) 
+2Ad 2 <%)u 2 (y A )] . 

Writing u 2 (y A ) as 

«2(z/a) = u 2 (yi,y 2 + 5(y)) + J d 2 u 2 (y^) d/J, 
and using the condition J 2 Xip(X) dX = yields 

|£>(«)(y)|<C [ 2 \D(u)\(y x )dX. 
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A careful analysis of the right-hand side then allows to prove that 

/ dm r iVl) dy 2 \D{u){y)\ 2 < C [ \D(u)\ 2 . 

JR Jmiu-T) JR Y 

For all additional details, we refer to |20j. 

Consequently, we have built an extension operator 

E : ^(Ry) i-> H x {R x (infw-r/,y)) 

such that for all u G H l {Ry) : 

\\ d ( u )\\lZ(r y ) < \\D(Eu)\\ L 2 {Rx{inlu; -, n)) < C\\D(u)\\ L 2 {Ry) , 

\\ u \\l 2 (R y ) - \\ Eu \\L 2 (Kx(mfw-r,)) < C\\u\\ L 2( RY y 

Second step: compactification of the problem. 

In the rest of the proof, we set Q := M x (inf uj — rj, Y). According to the first step, we now 
have to prove the existence of a constant C such that for any function u G H 1 ^) satisfying 
u ■ u\r = 0, 

\\u\\mQ) < C\\D{u)\\ L 2 {Q) . 
Of course, it is sufficient to prove that there exists a constant Ca such that for all k G Z, 

IMlL 2 (Q fe , A ) < C A \\D{u)\\ L 2 {Qk A) Vu G H l {Q kA ) s.t. u ■ v\ t a = (2.6) 

where 

Qk,A = QC\{y, kA< yi < {k + 1)A}. 

Assume by contradiction that ( |2.6| ) is false. Then there exists a sequence of relative 
integers (k n ) n >i and a sequence u n G ^(Qkn^i such that for all n, 

In the rest of the proof, we drop all sub- and superscripts A in Qk,A, Ik m or der to lighten 
the notation. 

Let v n := u n {- + (k n ,0))/\\u n \\ L 2. Then v n G ff 1 ((5o) for all n and 

v n ■ v\r kn = 0, KIU»(Q ) = 1> PK)||l 2 (Qo) ^ -■ 

According to the standard Korn inequality (see for instance [20]), there exists C > such 
that for all v G H 1 (Qq), 

\\Vv\\ L i m < C(\\v\\ L 2 {Qo) + ||D(v)|| L2(Qo) ). 

As a consequence, the sequence u n is bounded in ^(Qq). By Rellich compactness, there 
exists a subsequence (still denoted by v n ) and a limit function « G H 1 (Qq) such that 

v„ ->> v mw-H l (Q ), 
v n ^v mL 2 (Q ). 
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We deduce that D(v) = and \\v\\ L 2 
(C, i/o) £ (K x K 2 ) \ {0} such that 



1. Hence v is a non-zero solid vector field: there exists 



(2.7) 



v (y) = ( c y + Vo) for a.e. y E Q . 
On the other hand, for all n E N, for almost every y\ E [0, A], we have 

v n ,i(ik n (yi)h'k n ,2(yi) - vn^knivi)) 7L,i(yi) = °- 

Since the sequence 7fc n is bounded in W 1,co , up to the extraction of a further subsequence, 7fc n 
converges weakly - * in W l >°° towards a function 7. Since 7fc n ,i(yi) = yi for all n, we deduce 
that (7fe n — 7) • ei =0. We then pass to the limit in the identity (2.7) using the following 
facts: 



7fc„ — > 7 in L°°, and thus 

r-A 





)-t>»(7)l < C|l7fc„ -7lloo||Vw„| 



L 2 



0; 



• v n {pj) is bounded in i^ 1 / 2 ((0, ^4)), and thus 

Un(7)-»w(7) inL 2 (0,A). 

At the limit, we obtain 

(C7 + yo) • 7' = 0, 
|C7 + 2/0 1 2 = est. 



i.e. 



We deduce that 7 E 7*^, and thus 7*^ n {7^, k E Z} 7^ 0, which contradicts the assumption 
of the lemma. Thus (2.6) holds, which completes the proof. □ 



Remark 7. • PFe emphasize that condition (2.5) is probably not optimal. Indeed, (2.5) 
amounts to requiring that the inequality 

\\u\\ L 2<C\\D(u)\\ L2 

holds uniformly in each slice of length A. However, since our proof relies on compactness 
results in L 2 , it seems necessary to work in a fixed compact domain. Of course, if a 
more "constructive" proof were at hand (in the spirit of Lemma^), it is likely that (2.5) 
could be weakened. 



We have already pointed out that in the periodic case, conditions (2.5) and (2.1) are 
equivalent. In the general case, however, (2.5) is stronger than (2.1). Indeed, (2.1) 
merely requires the frontier T to be non-flat (uniformly on R), whereas (2.1) requires 
that is not invariant by rotation, in addition to being non-flat. 

We have used in the proof the following Korn inequality: since the function uj is Lipschitz 
continuous, there exists a constant Ck > such that 

IMItfi(fl) < Ck(\\u\\ L 2 {r) + \\D(u)\\ L 2 {R) ) \fu E H X (R). 

We refer to }2(tf (see also fllf) for a proof. The constant Ck depends only on the Lips- 
chitz constant of uj. The inequality holds without any assumption on the non- degeneracy 
of the boundary or on the behaviour of u at the boundary T. 
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We now prove Lemma [6j Assume that there exists Yq such that (Ky ) holds true. Let us 
first prove that (Ky) is true for all Y G (supw,Yo). Let u G H x {Ry) be arbitrary. Using a 
construction similar to the one of Nitsche (see [20]), we define an extension u\ G H of 
u such that 

Y\ = supw + 2(Y — supw), 
\\ui\\lz(r Yi ) < Ci\\u\\ L 2 {Ry) , \\D(ui)\\ L 2 {Ryi) < Ci\\D(u)\\ L 2 {Ry) . 

Iterating this process, we define sequences (Y n ) n >Q, (u n ) n >o such that u n G H l (Ry n ) and 
u n+ \ is an extension of u n for all n, and 

y„+i = supw + 2(Y n - supw), 

\\ u n+l\\L 2 (R Yn+1 ) < Cn+l|| u n||L 2 (fl y „), ||-D(«n+l) Hl^fly^) < Cn+1 ||-D(«n) IU 2 (fl r J • 

It can be easily checked that lin^^oo Y n = oo, and thus there exists no > such that Y„ > Yo- 
By construction, u no G H (Ry Q ) and there exists a constant C such that 

IKollL 2 (iiIy no ) < POOIIl 2 ^^) < C|l^( u )llL2(/? y )- 

Moreover, 

u = u no on i?y. 

> From (i£y ), we infer that 

IKoIIl 2 ^) < CVo ^(^o) IIl 2 (Ryo)' 

and thus 

IMIl^fly) < C\\D( U )\\ L 2 {Ry) . 

Hence (Ky) is satisfied. 

Let us now prove that (Ky) is also true for all Y > Yq. Let u G H l (Ry) arbitrary; then 
u G H l (Ry ), and 

IM| L 2 (i?Yb) < Cy \\D(u)\\ L 2 {RYo) . 

Moreover, according to the classical Korn inequality in the channel Ry , there exists a constant 
Ck such that 

h\\m(R YQ ) < C K {\\u\\ L 2 (Ryo) + \\D(u)\\ L 2 iRYo) J . 
Let S :=R x {Yq}. Then 

IMIl 2 (£) < c\\u\\ h i {Ryo) < C\\D(u)\\ L 2 {Ryo) . 
Now, for any y G R x (Y , Y), let y' G £ such that 

y' = y + t(l, —1) for some i G R. 

Then ^ 

u(y) = u(y') + t / (<9i - d 2 )u{y + *(1 - r)(l, -l))dr. 



Notice that 

(1, -1) • (d\ - 82) u = dim + 92^2 - (diU2 + &2U1), 
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and thus 



Similarly, 



\u{y)-(l-l)\ 2 dy < c(\\u\\ 2 +f j' Y \D{u)\ 2 
/R JY V JR Jy 

< C\\Du\\1 2{Ry) . 



°Ky)-(-l,-l)| 2 dy<C\\D( U )\\ 2 L2{RY) . 



Eventually, we obtain 

\Hl2(R y ) < C\\Du\\ 2 l2{Ry) , 

which completes the proof. 

3 Estimates for the no-slip condition 

In this section, we will prove Theorem [TJ In other words, we will establish, for any sequence 
A e £ [0,+oo], the well-posedness result and the error estimates that were established in [6 
for \ £ = 0. We will use the same general strategy, based on the work of Ladyzenskaya and 
Solonnikov [T7J. However, the handling of the slip type conditions will require new arguments, 
due to a loss of control on the skew-symmmetric part of the gradient. Moreover, we shall 
specify the dependence of the error terms with respect to <j). We shall of course put the stress 
on these new arguments. 

The starting point of the proof is an approximation scheme by solutions in truncated 
channels. Therefore, we introduce the notations 

Vk,l>0, U kji := UH{k < \xi\ < I}, U k := U_ k>k , 

for any set U of M?. We take as a new unknown 

v := u e — u°, u° := 1qu° 

where u° is the Poiseuille flow. As a new pressure, we take 

q := p + 12(pxi. 

It formally satisfies 

' -Av + u° -Vv + v- Vu° + v-Vv + Vq = l R e(-l2(f>,0), x G tt £ , 

div v = 0, x £ ST, 

v\ X2 =i = 0, / vi = 0, 

Jcr s 

v T \rt = A 6 {D{v)v) t \yz , v ■ v\rz = 0, 
where v is the inward pointing normal vector on T £ and 

V T = (Id - v ® v)v 



(3.1) 
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denotes the tangential part of v onT £ . The system is supplemented with the following jump 
conditions at the interface E: 

His = 0, [-D(y)e 2 + qe 2 ]\v = (-60,0). 

In order to build and estimate the field v, we consider the approximate problems in Q, e n 

' -Av + u° -Vv + v- Vu° + v ■ Vv + Vq = 1^(120,0), x E U £ n , 



\xi=n 



\xi=—n — 

v T \rz 



div v = 0, x G f^, 
v\ X2 =i = 0, 

= \ £ (D(v)v) T \ r e v ■ u\ T e 



(3.2) 



and 



N| En = 0, [-D{v)e 2 + qe 2 ]\v n = (-60,0). 
The proof divides into four steps: 



1. We construct a solution v n of the approximate system (3.2). 

2. We derive H^ loc estimates on v n that are uniform in n. This yields compactness of 
(v n ) n , hence, as n goes to infinity, a solution u of (NS e ). 

3. We prove uniqueness of this solution u. 

4. We deduce from the previous steps the desired 0(y/e) bound in Hh for u — u°. From 



there, using duality arguments, we get the O(e) bound in L 



uloc' 



Step 1. The wellposedness of (3.2) relies on an a priori estimate over Q e n . Multiplying 
formally by v, we obtain 



/ \D(v)\ 2 + {X s )- 1 I' \v T \ 2 = - [ (u° <g> v + v ® u°) : D(v) - [ 124> Vl +[ 



where v T denotes the tangential part of v. As HVn ^ < C 4>, we obtain 

\\ D t v )\\h{n%) ^ c $ {\\ v \\L^n^\\D{v)\\ L 2 (nik) + ^/n£\\v\\ L 2 (R e n) + V" IMIza(£«)) ■ (3-3) 

As v is zero at the upper boundary of the channel, Poincare inequality applies, to provide 

\Hl^) < CUVull^n.). 

where C depends only on the height of the channel. Let 

n e bl := {x, x 2 > euj(xi/e)} 

the "rough" half plane, and v £ H (flf,) the extension of v which is zero outside f^. We can 
apply to v the results of Nitsche |20] on the Korn inequality in a half plane bounded by a 
Lipschitz curve: one has 

||V?)|| L 2 ( ^ ;) < C\\D{v)\\ L 2 {nii) 
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where the constant C only depends on the Lipschitz constant of the curve. For fif,, this 
Lipschitz constant, and therefore the estimates, are uniform in e. We insist that this inequal- 
ity is homogenenous: it does not involve the L 2 norm of v, contrary to the more general 
inhomogeneous Korn inequality. Back to v, we get 

UVvll^n.) < C\\D(v)\\ LHn s n) . 

Hence, denoting for all k £ N 

E k ■= \Hh(n%) + W^ v \\h{ni) + \\ D { v )\\h{ai) 
the combination of Poincare and Korn inequalities leads to E n < C\\D(v )\\ 2 ^2/^ey 

Now, by rescaling either the Poincare inequality when A e = 0, or the inequality in (H') 
when X £ ^ 0, we get 

IMIz^) < Ce\\Vv\\ L 2 {R s n) . 

Then, we deduce 



Back to the energy estimate (3.3), we end up with 



E n < C \\D(v)\\l 2(nEn) < CXH^e) \\D(v)\\ L 2 m + V^\\Vv\\ L 2 (m ) 

, \ I (J'2 

< C'(f>[E n + ^e^EA < C'(f)E n + -E n + — <\? ne. 



2 2 

so that, for <f) small enough, we have the global estimate 

E n < C(p 2 ne. (3.4) 

Thanks to this estimate, one obtains by classical arguments a variational solution v n £ 



H l {Q, e n ) of (3.2). The uniqueness of this solution (for <ft small enough) is deduced from the 
same kind of energy estimates, performed on the difference of two solutions. We leave the 
details to the reader. 

Step 2. The next step in the proof of Theorem [l] is the derivation of uniform H^ loc bounds 
on v n - The idea, which originates in a work of Ladyzenskaya and Solonnikov, is to prove by 
induction on k' = n — k that 

E k < C (j) 2 (k + 1) e, C large enough. (3.5) 

Once the bound on the E^s is proved, we can use it with k = n — 1, so that E\ < C (j) 2 e. 
This gives a control on a unit slice of the channel around x\ = 0. But as will be clear from 
the proof of the induction relation, x\ = plays no special role: in other words the same 
bound holds for any unit slice of the channel, which gives the uniform Hh bound. 



Let us now describe the induction process. First, by (3.4), the induction assumption holds 
with k' = 0. To go from k' — 1 to k' , that is from k + 1 to k, we shall need the following 
inequality: 

Vk < n, E k < Ci ( E k+1 - E k + (E k+l - E k fl 2 + <p 2 (k + 1) s) . (3.6) 
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This inequality at hand, and assuming E k +i < Co 4> 2 (k + 2) e, one obtains straightforwardly 
that Ey. < Cq 4> 2 (k + 1) e for all k > Co — 1 provided Co is chosen large enough. For 
k < Co — 1, we have merely < Co (j) 2 {[Co\ + 1) e . Hence, up to a new definition of the 
constant Cq, we obtain inequality (|3.5|). 



An inequality similar to (3.6) has been established by one of the authors in [6], for the 
case X s = 0. The discrete variable k is replaced in [B] by a continuous variable r/, but the 
correspondence from one to another is obvious. We also refer to the original paper [IT], and 
to the boundary layer analysis in [13], in which similar inequalities are derived. 



Relation (3.6) follows from localized energy estimates. We introduce some truncation 



function Xk = Xfc( x i)> such that Xk = 1 over Sl|, Xk = outside and |x' fc | < 2. 

Multiplying by Xk v within (3.2) and integrating by parts, we deduce that 



/ Xk\D(v)\ 2 + (XT 1 [ Xk\Vr\ 



< 



Xk(u° ®v + v®u°) : D(v) 



12cpXkVi + / 6<MXfc 

D(v) : (Vxk ® v) + / {u° v + v (E> u°) : (Vxk ® «) 
<» ./<> 

7 



iff iff - =1 

where V^fc = (Xfc>0)- The r.h.s. of the inequality has two different parts: 

• The first three terms are very similar to those of step 1. They are treated along the 
same lines: 

3 

< C4> (E k+l + y/{k + 1) e y/E^l 

3=1 



The remaining terms involve derivatives of Xk- they are supported in Slk,k+i- Standard 
manipulations yield the bounds: 



1-^4 1 + 1-^5 1 < C(Ek+\ — Ek), \h\ < C(Ek+i—Ek) 



3/2 



• The treatment of the pressure term is a little more tricky. We decompose 



i %k+i u Si fc,fc+i> 



n k%i ■■= n|,fc+in{±xi>o}. 



The zero flux condition on u implies that f e,± f(xi)vi = for any function / depend- 
ing only on x\. Thus, 



h 



' Xk v l 



k,k+l 



(q-q k )x k vi + 



k,k+l 



(q-qt)x'k v i, 



h,h+l 



where is the average of q over fi^jfe , i" Then, we use the well-known estimate 



L2(0) 



< C\\Av + f\\ H -x 



(©) 
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for the Stokes system 

-Av + Vq = f, div u = in O 

where C only depends on the measure of O and the Lipschitz constant of dO. We take 
here O = fi^, j (so that the constant is uniform in k and e), and 

/ : = -div (u° ®v + v®u° + v®v) +1^(120,0). 
From there, one gets after a few computations 

\Ir\ < ft ll^nf^ + I'*" ^'^(nt.u) ll w M<W) 



< C \ 4> \/e y 7 E k+ i — E k + {E k+ i — E k ) + (E k+ i — E k 

< C U 2 e + (Sfc+a - £ fc ) + (E k+1 - E k f' 2 



We refer to [6] for more details. By gathering all the inequalitites on the iVs, we obtain for 
cj) small enough 



J Xk\D(v)\ 2 < C (E k+l - E k + (£ fc+1 - E k f' 2 + (ft 2 (k + l)e 



+Cct> [E k + ^{k + l)e^E~ k 



(3.7) 
(3.8) 



Now, we have 



I / |2| 2| 



Xk\D(v)\' > / x^l^WI 2 > / \D( Xk v)\ 2 
> f \D( Xk v)\ 2 -4(E k+1 -E k ). 
As XkV is zero outside 0|, 1; we can proceed as in Step 1, to get 

f |^(Xfc^)| 2 > cE k+1 > cE k + c(E k+1 -E k ). 
J ri e 

Finally, 

f Xk\D(v)\ 2 > CE k - c(E k+1 -E k ). 
iff 

Combining this inequality with (|3.7|) gives the result provided 4> is small enough. 



As we have already explained, once inequality (3.6) is proved, one obtains easily a 0(y/e) 
Huioc bound on v n . By standard compactness arguments, any accumulation point v of (v n ) 



is a solution of (3.1). It provides a solution u £ of the original system (NS e ). Moreover 



\U 6 — It !! TTl 



< C(ft^~E, 



(3.9) 



iSiep 5. It remains to prove the uniqueness of the solution in H^loc- 

Let now v = u £ ' 2 — u 6 ' 1 the difference between two solutions of (NS e ). It satisfies 



-Av + div (u 6 ' 1 ®v + v® if' 1 + v ® v) + Vq = 
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together with div v = and homogeneous jump and boundary conditions. We can always 
assume that u 5 ' 1 satisfies the bounds in (3.9). Then, performing energy estimates similar to 
those of Step 2, we get, for <f> small enough: 

Vfc, E k < C ({E k+1 - E k ) + (E k+1 - E k f' 2 + 4> ^~sE k 



We have used implicitly that 



XkU 6 ' 1 ® v : V(u) 



< 



eA \ M \V(v) 



3=0 Jn j,j+i 



iJ+1 



k 

< C(/>^E k+1 . 

3=0 



^ H u£,1 |l^ oc I^ll u llHi(^, i+ i) 



As v belongs to H^ loc , E k+ i — E k is bounded uniformly in k: eventually, for <\>yfe small 
enough, we get E k < C for all k, which means that v is of finite energy. The fact that v = 
then follows from a classical global energy estimate, performed on the whole channel Sl e . This 
concludes the proof. 



Step 4- Note that, by the previous steps, we have established not only the well-posedness, 

H 1 , (Q) 



but the H^ loc estimate 



\u £ - u°\\m m\ < C 4>y/e. 



From there, one obtains that 



\u £ -u°\\ L 2 (s) < C4>e. 



The L^ oc (f2) estimate follows from estimates on a linear problem in the channel O: 

-Av + Vq + v° ■ Vv + v ■ Vu° = div F £ , 

div v = 0, 



where v = u £ — u°, ip £ = v\% = 0(<pe) in E 2 uloc and F £ = v®v = 0{4>e) in L^ loc thanks to 
the H^ loc bound. By a duality argument, one can then prove that ||v||^2 ; ^ is also 0(4>e) in 
E 2 uloc . This duality argument is explained in the paper (BJ section 3.2]: as the slip condition 
in the boundary condition at d£l £ plays no role, we skip the proof. 

4 Boundary layer analysis 

In order to improve our description of u £ , we must analyze the behaviour of the fluid in the 
boundary layer. The starting point of this analysis is a formal expansion: we anticipate that, 
near the rough boundary, we have 

u £ {x) = u°(x) + 6<pev(x/e) 
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where u° is the Poiseuille flow, and v = v(y) is a boundary layer corrector, due to the fact that 
u° does not satisfy either the Dirichlet boundary condition when A = 0, or the slip boundary 
condition when A 7^ 0. We shall focus on the latter case, which is the new one. Classically, the 
rescaled variable y belongs to the bumped half plane Q, bl := {y, y% > co(yi)}, and by plugging 



the expansion in ( NS e ) , one finds that 



(D(v)v) 



-Av + Vp = 0, 
div v = 0, 
= -(D((y 2 ,0))u) T , 
v . v = -(y2,0) • V, 



yen bl , 
y e n bl , 
y e an bl , 
y e an bl . 



(BL) 



where 



v(y) :-- 



1 



= (-^(lft),l) 



Vi + 7 /2 (yi) 

is a unit normal vector. The inhomogeneous boundary terms come from the Poiseuille flow 
(^2(1 — X2) ~ ey2 near the boundary). 



System ( BL ) is different from the boundary layer system met in the former studies on wall 



laws: it has inhomogeneous Navier boundary conditions, instead of Dirichlet ones. Neverthe- 
less, we are able to obtain similar results, as regards well-posedness and qualitative issues. 



4.1 Well-posedness of the boundary layer 

First, we have the following well-posedness result: 



Theorem 8. System (BL) has a unique solution v G Hi oc (Q, bl ) satisfying: 



sup 
k Jn bl 



\Vv\ 2 < 



-00 



where for all k,l, 



tt b ir\{k < yi < I}. 



k,k + l 



The proof follows closely the lines of [13J, where the case of an inhomogenenous Dirichlet 
condition (instead of Navier) was considered. The main difficulty comes from the unbound- 
edness of the domain, which prevents from using the Poincare inequality or assumptions like 
(H') or (H"). To overcome this difficulty, there are two main steps: 



1. One replaces system (BL) by an equivalent system, set in the channel 



n bl n{ y2 <o}. 



This equivalent system involves a nonlocal boundary condition at y 2 
Dirichlet-to-Neumann type operator. 



0, with a 



2. Once brought back to the channel Q b1, , one can follow the same general strategy as in 
the previous section, based on the truncated energies 



+ \\m\\^-y 
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Let us give a few hints on these two steps. 

Step 1. It relies on the notion of transparent boundary conditions in numerical analysis. 



The formal idea is the following: the solution v of (BL) satisfies the boundary value problem 

(4.1) 



-Av + Vq = 0, y 2 > 0, 
V • v = 0, y 2 > 0, 

V \y 2 =Q = VQ, 



where vq = v\ V2= q. Using the Poisson kernel for the Stokes problem in a half-plane, we have 
the representation formula: 



v(y) 



G(t,y 2 )v (y 1 -t)dt, q(y) 



Vg(t,y 2 ) -voim -t)dt 



(4.2) 



where G(y) 



2y 2 



v(yl + y 2 2 ) 2 yl ' ' 9[y> 



2y 2 



K{yl + y 2 2 Y 

Thanks to this representation formula, we can express the stress 

(2D(v)u - qv)\y 2 =o = -2d 2 v + (d 2 vi - <9if 2 )ei + qe 2 
in terms of v at y 2 = 0. Formally, this leads to some relation 

(-2D(v)e 2 + qe 2 )\ {y2=0} = DN(v\ {y2=0} ) 
for some Dirichlet-to-Neumann type operator DN. Hence, and still at a formal level, we can 



replace the system (BL) by the following system in Q b 



-Av + Vq = 0, ye IT", 

v-u = o, yen bl -, 

(D(v)u) T = -(D((y 2 , 0))v) T , y € dn bl '~ , 
v ■ v = -(y 2 ,0) ■ u, yedn bl '-, 
{-2D(v)e 2 + q e 2 )|{ y2 = } = DN(v\ {y2=0} ). 

A rigorous version of these formal arguments is contained in the next proposition 



(4.3) 



Proposition 9. (Equivalent formulation of (BL)) 
i) (Stokes problem in a half-plane) 

For all vq 6 H^ 2 c (R) there exists a unique solution v G /Z^ C (M^_) of (4.1) satisfying 

rk+l r+oo 
/() 



sup 

fcez Jk 



\Vv\ dy 2 dy\ < +oo. 



(4.4) 



ii) (Dirichlet-to-Neumann operator) 
There is a unique operator 



DN : Hl /2 C (R) ^ V'(R) 
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that satisfies, for all Vq G H^JR), and all ip G C^°(R\) with V • ip = 0, 



D(v)-D(cp) =<DN(v ), P\ {y2=0} >. 



(4.5) 



where v is the solution of (4.1). Moreover, for all vo 6 H uloc (R), the operator DN(vq) 
can be extended to a continuous linear form over the space Hl /2 (R) of H 1 / 2 functions 
with compact support. 

iii) (Transparent boundary condition) 

Let (v,q) be a solution of ( |BL[ ) in Hi oc (Q bl ) with sup fc f n m \Vv\ 2 < +oo. Then, it 



satisfies (4.3). 

Conversely, let v~ in H}., (Q ~) be a solution of (4.3). Then, the field v defined by 



v~ inn bl ~, v := I G(yi - t,y 2 ) v~ (i,0) dt for y 2 > 
Jr 

is a solution of (BL) in Hl oc {Q bl ) such that sup fc j^bi |Vw| 2 < +oo. 



k,k + l 



Proof of the proposition. The proof of the proposition is almost contained in [13]. The 
only difference lies in the definition of the Dirichlet-to- Neumann operator. In |13j . the full 
gradient is used in the definition of DN, instead of its symmetric part. Here, in order to adapt 
to the Navier condition at the rough boundary, D(u)v substitutes to d u u, and, subsequently, 
(4.5) substitutes to the relation 

f VvVp =<DN(v ), ip\ {y2=0 }> ■ 

used in [13 1 . As these minor changes do not play any serious role, we skip the proof. 

Step 2. By the previous proposition, in order to prove well-posedness of the boundary 
layer system, we can work with the equivalent system (4.3 ). As it is set in a bounded channel, 
it is amenable to the kind of the analysis performed in the previous section, for the study 
of the no-slip condition. The keypoint is again to have an induction relation between the 
truncated energies. However, the nonlocal DN operator prevents us from deriving a local 
relation like (3.6). We are able to show the following more complicated relation: there exists 
r] > such that, for any m > 1, 



k + 1 H sup (Ej + \ — Ej) + m sup {Ej+i 

TTl^ jy} t j rm k+m>j>k 



E 3 ) 



(4.6) 



The same relation was established in |13j . for the boundary layer system with a Dirichlet 
condition. The proof starts with the the change of unknowns u := v + (y2,0), p := q, that 
turns (4.3) into 



-Ait + Vp = 0, 
V -u = 0, 
{D{u)u) T = 0, 
u ■ v = 0, 



yen bl >~, 
yen bl >-, 
y e dn bh 
y e OQbi, 



(-2D(u)e 2 + qe 2 )\ {y2= Q } = DN(u\ {y2=0} ) - (1,0). 
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Afterwards, energy estimates are performed, testing against XkU. The only change with 
respect to [13] due the Navier condition is the treatment of the lower order terms. When a 
Dirichlet condition holds at the boundary, one can rely on the Poincare inequality, to obtain 

E k < C If u _ |Vu| 2 + E k+l -E k \ 

and then to control E k from the energy estimate (which gives a bound on the gradient only). 
This is no longer possible in the case of the Navier condition. Moreover, we cannot proceed as 
in the previous section, using the Dirichlet condition at the upper boundary of the channel. 
Indeed, in our boundary layer context, a non-local condition holds at the upper boundary. 
This is where the assumption (H") is needed: it easily implies that 

E k < C (J u \D{u)\ 2 + E k+X -E k \ 

and can therefore be controlled from the energy estimate (which gives a bound on the sym- 
metric part of the gradient only). From there, all computations and arguments are similar to 
those of |13j . We refer to this paper for all necessary details. 



4.2 Qualitative behaviour at infinity 



As the solution of (BL) is now at hand, we still need to show its convergence to a constant field 
as y2 goes to infinity. Here, some ergodicity property condition must be added. We consider 
the stationary random setting: we take uj to be an ergodic stationary random process (on a 
probability space (M, /i)), obeying the assumptions of Theorem[2] We then state the following 
proposition: 



Proposition 10. There exists a £ M. such that the solution v of (BL) satisfies 

v(y) ->■ (a,0), as y 2 ->■ +oo, 
locally uniformly in y\, almost surely and in L P (M) for all finite p. 



This proposition is based on the integral representation (4.2), and the ergodic theorem. 



It has been proved in [6j (the condition at the rough boundary does not play any role). 

Remark 11. It is possible to derive upper and lower bounds for the "slip length" a. Such 
bounds were established by Achdou et al. in the periodic setting. They are still valid in the 
random stationary case. Along the lines of J3\j, one can prove: for uj £ C 2 (M) n VF 2 '°°(IR), 

Y m ax ^ Q! ^ Y mm , Y m i n . — inf UJ, Y max . — supcj. 

K R 

In order not to burden the paper, we skip the proof. 
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5 Estimates for the Navier condition 



This section is devoted to the proof of Theorem [2} The main novelty lies in the derivation 
of almost sure estimates. Indeed, to our knowledge, the previous convergence results dealing 
with a stationary ergodic setting were all stated in a norm involving an expectation (see 
j6j H2]). The main steps of the proof are the same as in |13} [6]: the idea is to build an 
approximate solution which consists of the main term u , the boundary layer corrector v, and 
two additional correctors u 1 and r 6 . We will review briefly the definition and well-posedness 
of u 1 and r £ , and focus on the estimates which are required for the proof of Theorem [2| 

• We start with some regularity estimates for the function v which solves the boundary 
layer problem (BL): 



Lemma 12. Let j3 G N be arbitrary, and let v be the solution of (BL). Then for all a > 0, 
there exists a constant C , depending only on the Lipschitz constant ofu, on (3 and on a, such 
that 

rk-\-\ roc 2 

sup / / V p Vv < C. 

k& Jk J a 

In particular, v G L°°(M x (a, oo)) for all a > 0.. 

Proof. The arguments are the same as in |13l Proposition 6]. According to Propositi on |9| 
Vv G L^ loc (Q bl ), and vq = G L^ oc (M). Since v is given by the representation formula (4.2) 
in the upper-half plane, by differentiating under the integral sign in (4.2), we obtain 

fk+l 



V^Vv 



< 



< 



< 



k+l 



k+l 



V l3 VvG(t,y 2 ) \v {yi 
1 



dt 



t 2 + y 2 2 



vo(yi -t)\dt 
k+1 r Myi- 



t 2 + y 2 2l '" i;/l 



t)\ dt 

2 



dyidy 2 



dyidy 2 
t)\ 2 



t 2 + y 2 2 



dt dyi dy 2 



< Cp\\v Q \ 



1 



1 



2/2 i 2 



-^dtdy 2 < C/3. a \\vo\\ L 2 , 

y* ulocV 



□ 

• We now prove the following result, which is crucial with regards to the derivation of 
almost sure estimates, and which is the main novelty of this section: 



Proposition 13. Letv be the solution of the boundary layer system (BL). Then the following 
estimates hold almost surely as e — > 0: 



v 2 



x[ 1 

J — 

e e 



K7 £ ) - ( a ,°)\\L?(n R ) 



dx'i 



+ 



H 3 (-R,R) 



Vl 



x± 1 

£ £ 



a 



H 3 (-R,R) 



o(l). 
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Remark 14. This result combines two main ingredients: the deterministic construction of 
the preceding section, which eventually led to Lemma [M and the almost sure convergence 



of the corrector v in the stationary ergodic setting (see Proposition^!^. We emphasize that 
both items are important here. In particular, it does not seem possible to prove almost sure 
estimates by using a probabilistic construction of the boundary layer as in f^- 

Proof. We use an idea developed by Souganidis (see |21j). Let 5 > be arbitrary. Then, 
according to Egorov's Theorem, there exists a measurable set Ms C M and a number ys > 
such that 

\v(0,y 2 ,m) - (a, 0)| < 8 Vm£ M s , Vy 2 > ys, 

P(M!) < 5. 

Without loss of generality, we assume that ys > 1. 

Now, according to Birkhoff's ergodic Theorem, for almost every m there exists ks > 
such that if k > ks, 

A s = A s {m) := {y\ G M.,T yi m € Ms} satisfies: \A S n (-k, k)\ > 2k(l - 25). 
For all R > 1, e > 0, we have 

R-l\\v(./s)-(a,0)f L z(n R ) = i [ R/e [ 1/E \v(y,m)-(a,0)\ 2 dy 



R 

2 rR/e rl/e 

/ / U,((\ <ll„ n- m\ — (r* r^^ 1 2 

R 

c 2 rR/e ry s 

/ / u,rn o, n t m\ — (™ flV 2 

R 



R/e Jw(yi) 

\v(0,y2,T yi m) - (a,0)\ z dy 



R/e Ju(jn) 

\v(0,y 2 ,T yi m) - (a, 0)^ dy 

R/e Ju(yi) 



_2 r R/e rl/e 

£ I I , ,„ s , -,|2 



-R/e Jys 



+ ~R ' i lr vi m ^ M s l u (°>2/2,r yi m) - (a,0)\ dy 



£ 2 rR/e rl/e 

+ / ir^meM? \v(0,y 2 ,T yi m) - (a, 0)| dy 

K J -R/e Jys 



3 



3=1 



We have clearly, by definition of As, 



1 2 

h < esu > P^7 Mv) ~ ( a M\tf{((-R<,Ri)x(-l,y s ))r>nU) < °s 



, I | i; t I — - ' 

/ 2 <^ 2fil 5 2 <2<5 2 . 



Notice that the constant Cs in the first inequality depends on the random parameter m. 

As for the third integral, recall that v S L°°(]R x (1, oo)) according to Lemma 12 and that 
ys > 1- Thus we have, for all R > 1 and if e < 1/ks, 

h < 45 (\a\ 2 + lid' 2 



IL°°(Kx(1,do)) 

Gathering the three terms, we deduce that the first estimate of the Lemma holds true. 



2N 



Concerning the second estimate, we define the vector field u[y) = v(y) + (y 2 ,ti). Notice 
that u is divergence free and that u ■ u = at the lower boundary of fl bl . Consequently, 
following [6j, Proposition 14, we write 



dx 1 = e I u 2 [ yi, - ) dyi 



ui(—,y2)dy 2 -s u 1 (0,y 2 )dy 2 



l/e 
u(xi/e) 

(vi — a) ( — ,2/2) dy 2 — s / (vi - a) (0, y 2 ) dy 2 (5.1) 

u(xi/e) - 1 ■/'•''•nl 

^ 2 (^)-- 2 (o) 



+ a [oj 



l/e 



The last term is bounded in L°° by edMI^ + SUcjUoo), and thus converges towards zero in 
the appropriate norm. As for the other two terms, set 

1/6 



U e (yi) :=e (vi - a)(yi,y 2 )dy 2 . 
Using the same decomposition as previously, we write, for 5 > arbitrary, 



U £ (yi) 
and thus 



rys rl/e /-l/e 

/ (vi - a) + £l Tyi ueA s / («l - a) + elT S1 we^ / («i - a) 



!),) 



U £ (yi)\ <e \vi-a\ + 8+ \\vi - a\\ Lac(Rx(1:0o)) l TyiUje A c s - 



Consequently, we obtain, for all i? > 1 and for e small enough (depending on 5), 



1 



/ t/ e ( — ) d Xl < C~(y d -infu) / / 
J-R v e 7 J-ii/e JuOyi 



\v\ — a\ 



-R/e Jw{ Vl ) 

+C (S 2 + 48\\vx - a||£oQ(Rx(i,oo)) 
< C s e 2 + C5. 



Hence, as e vanishes, 



1 



sup- 

R>1 tt J-R 



R 



U £ (^r)\ dx x = o{\). 



The second term in (|5.1|) is easily treated: since it does not depend on x-y, we have 

i-R r l/e 



1 

SU P ~F> , 

R>1 J-R 



w(0) 



(«1 - a) (0, 2/2) dys 



l/e 



(«i - a) (0,2/ 2 ) dy 5 



w(0) 
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Since (vi — a) (0,2/2) vanishes almost surely as y 2 — > 00, we infer that 

•l/e 



lim e / (yi — a) (0, 2/2) ^2/2 = a.s. 
-Mo) 



This proves that 



as e vanishes. 
The estimates 







L 2 (-R,R) 



o(l) 



?'2 



1 

SU P m72 



Xl 1 

e e 

X\ 1 

e e 



L 2 (-R,R) 



a 



oil) 



L 2 (-R,R) 

are derived in a similar fashion. There remains to prove that for 1 < k < 3 



,-k 



X1 1 



0(1). 



L 2 (-R,R) 



Notice that it suffices to prove that for k = 1,2,3, 



lim yWv(0 , 2/2, m) = almost surely for m E M. 



(5.2) 



Then the same arguments as above allow us to conclude. 



In order to obtain (5.2), we use the same estimates as in [6], Proposition 13. We write 

'l M 



Hi 



;(0,ite,m) = Jj l d*+ 1 G(-y[,y 2 ){^J V \v Q {-z 1 m) - (a,0))dz\ dy[. 



Since G is homogeneous of degree —1, it can be easily proved that for k > 1, 



d k Vl G( yi ,y 2 ) <C k (yl + y 2 2 ) 



yid^ 1 G(y 1 ,y 2 ) <C k (yi + y z 2 



2\- 



Now, let 5 > be arbitrary. Almost surely, there exists ys > (depending on m) such that 

1 f yi 



2/1 Jo 



(vo(z,m) — (a,0))dz 



<$ if I2/1I > 2/<5- 



As a consequence, 

/ 2/i^' 1 +1 G(-y 1 ,2/ 2 
< C k 5 / (2/i 2 + 2/|) _ ^ i rf2/i 



1 pi 



2/1 Jo 



(v (-z,uj) - (a,0))dz dy[ 



< C k 



s_ 
2/2 
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On the other hand, 



J\yi\<ys \ViJo 



< c k 



Vi Jo 

\(vo(-z,m) - (a,0))| dz 



o 



(y (-z,m) - (a,0))dzj dy[ 

{yl + ylT^dyi 

\yi\<ys 



< C 5 



it 1 ' 



Gathering the two terms, we infer that for all 5 > 0, there exists C$ > such that 



y 2 d^v(0,y 2 ,m) 



<5 + 



Cs 

2/2 



Vy 2 > 0, 



and thus the quantity in the left-hand side vanishes almost surely as y 2 — > oo. 



□ 



We are now ready to prove the convergence result stated in Theorem [2j Following |13j . 



we set 



u lp P ( x ) '■= u °( x ) + 6^ (j) + ev^ix) + er e {x) + 60eV (^j , 



where the correctors u 1 and r e ensure that u| pp satisfies the Dirichlet boundary condition at 
the upper boundary and the zero flux condition. The term v 1 is a boundary layer term which 
compensates the tangential trace of u° + 6(j)ev at the rough boundarjQ Additionnally, u 1 is 
intended to be 0(1) while r £ = o(l). 
> We choose u 1 to be the solution of 

-Au 1 + u° ■ Vu 1 + u 1 ■ Vu° + Vp 1 = 0, if!!, 
V • u 1 = 0, 



u 



|x 2 =0 



0, u\ 



1*2 = 1 



-6(f>a, 0), 



on . 



Notice that we assume that u 1 satisfies a no-slip condition at the lower boundary. This stems 
from the non-degeneracy of the frontier T £ : in order that the non-penetration condition is 
satisfied at order e, u 1 must vanish at x 2 = 0. We recall that the same argument led to the 
no-slip condition for u° at x 2 = 0. 

Hence the vector field u 1 is exactly the same as in [13] and is a combination of Couette 
and Poiseuille flows: 

u\(x) = 6(f)(— 4ax 2 + 3ax2), u\(x) = 0, x G 

and we extend by zero outside Q. 

> The additional boundary term v 1 solves the system 

-Av 1 + Vqi = in n bl 
v l • v = (2/2 >°) • v on r ' 

A°(D(u 1 )i/) r = (u + (y 2 , 0)) T + A°( J D((y 2 2 , 0))^) T on I\ 



1 It can be checked that this extra boundary layer term is needed only when the original slip length A e is 
such that A < 1. 
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under the condition 



In 1 , 



sup 

k&i Jag.,, 

fc,fc+l 



\Vv 



1 1 2 



Using the same techniques as in Sectionpl energy estimates in H^ loc (U_) for v 1 can be proved, 
leading to existence and uniqueness of it. Additionally, there exists /3sM such that 

lim v 1 (y 1 ,y 2 ) = (/3,0) 

1/2-5- OO 

almost surely. 

> As for r e , we use the following Lemma: 

Lemma 15. There exists a vector field r e £ Hf 0C (Cl) satisfying 
V • r e = 0, 

r k=o = 0, rf X2=1 = 60 [(a, 0) - v (f , I) - ev' (f , 1)] , 
J CT r f = " J^v u i - 6 <^ + ev\)(x/e) + 60a, 



and such that 



1 

R>1 R 1/z 



l^ 2 (^-fl,fl) = and Il r lvy2>°°(n) = 0(1)- 



The Lemma follows directly from Proposition 13 and the construction in Proposition 5.1 
in [6]. Once again, we extend r e by zero outside f2. 
By construction, the function u| pp satisfies 

f -A< pp + < pp • V< pp + V/ app = div G £ + f in ft e \ S, 



app 1 ""app 
ul Lpp\x2=l = 0; 

< PP • ^lr- = 0, 



«pp)rb = A°( £) «pp) I/ )r|r- +9 £ , 

s K PP ]|s = 0, [D{ul pp )e 2 -p £ e 2 ] |s = D(er £ + eu 1 )e 2 |s =: ¥> e 



where 



G £ = eu u ® (^60 ( v 
+ e(6<j>(v 



(a,0) + SV 1 



+ r c 



(q,0) +r £ j 



+ e z ( 60 ( v ( -J + eu 1 yj + r £ J ® (^60 (v ^ 



£«* I - ) ) + r 



and 



£ V(x/ e )-(x^,0)) T | 



According to the estimates of Section [4} Proposition |13| and Lemma 15 we have 

1 



1 



lk e ||i2 (Sjl ) = o(e) + O(e0), 
sup|b £ || L 2 (r| ) = 0(e 2 ). 
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Consequently, setting w £ = u e — uf DP , we obtain 



Aw £ + (u £ ■ V)w £ + (w £ ■ V)< pp + Vq £ = -div G £ - f inft £ \ S, 



(5.3) 



and u; £ satisfies the same boundary and jump conditions as u^ pp . 

The next step is to derive energy estimates for the above system. The proof goes along 
the same lines as the one in Section [3j and therefore, we skip the details. The main steps are 
the following: 



1. First, we derive an energy estimate in f2 £ for a sequence (u> £ ) rae ^ satisfying (5.3) in VL £ n 
with homogeneous Dirichlet boundary conditions at x\ = ±n; more precisely, we prove 
that for cj) small enough, 



no(e 



[ (|l»«)| 2 + ki 2 + |v< 

2. By induction on k, we prove that for all n > 1, k G {1, • • • , n — 1}, 

/ (| J D«)| 2 + K| 2 + |V<| 2 ) = A; ( e 2 ). 

3. Passing to the limit as n — > oo, we deduce that for all k > 1, 

(\D(w £ )\ 2 + \w £ \ 2 + \Vw £ \ 2 ) = ko{e 2 ). 



There are two main differences with the estimates of Section [3l The first one lies in terms 
of the type 

|(«-VK pp ) -<|; 

indeed, because of the boundary layer term v, Vm| pp does not belong to L°°(Q £ ) in general. 
Therefore, using Sobolev embeddings, we have 



Vv 



fc-l 

dx = J2 

j=—k " "jj+i 
k-i 

j=-fc 



K(aO| : 



e |4 



1/2 



dx 
i+i -i/ e 



1/2 



|V^(y)| 2 o!y 



e i|2 



The second difference comes from the boundary term g £ , namely 

J re Xk(w £ n ) T g £ <j^J re Xk\{w £ n ) T \ 2 + C{k + 1) " 



£ 



The first term of the right-hand side can be absorbed in the boundary term coming from the 
integration by parts of J Au; £ ■ w £ n . The second one is clearly o((k + l)e 2 ). Notice that this is 
the reason why we need the additional boundary layer term v : if we merely take 



u 



app 



u +ev [ — ] + eu (x) + er (x) 
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then g £ = 0(e), and the second term in the right-hand side of the preceding inequality is 
O((A; + l) e 2 /A ). 

The inequality relating E}. and E k+ \ (with the same notation as in Section [3| becomes in 
the present case 

E k < cf^Ek^-Ek + ^Ek^+ev^Vk + lyfE^+ik + iy 
+ C(0 + e V (e)Vk + l)(E k=1 - E k f' 2 , 

for some function 7] such that lim + rj = 0. By induction we infer easily that 

E k < ke 2 rn(e), 

for some other function rji vanishing at zero, which completes the second step described above. 
The two other steps are left to the reader. 
We infer that 

sup — jy^ 1 1 u £ - < pp \\ Hl{n e R) = o(s) almost surely. 



On the other hand, let u N be the solution of (NS )-(Na) with A = 6<pae. Then the function 
u N is explicit: as in [13], we have 



u N = (fyU N (x 2 ),0) withU N (x 2 ) = -^^x 2 2 + — i —x 2 + £a 



1 + \ea 1 l + Aea 1 + 4ea ' 
so that 



N _ „,0 



u° + 6<pe(a, 0) + eu l + 0(e 2 ) in L 2 uloc (tt). 



From there, we obtain 



\\ uN ~ n a PP || jL 2 ( ^| i) = o(e) almost surely. 



Theorem [2] follows. 
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